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£Nj . In this series of three papers, we generalize the derivation of dual photons and monopoles by 

Polyakov, and Banks, Myerson and Kogut, to obtain approximative models of SU(2) lattice gauge 
theory. Our approach is based on stationary phase approximations. 

In this third article, we start from the plaquette representation of 3-dimensional SU(2) lattice 
gauge theory. By extending a work of Borisenko, Voloshin and Faber, we transform the expectation 
value of a Wilson loop into a path integral over a dual gluon field and monopole variables. The 
action contains the tree-level Coulomb interaction and a nonlinear coupling between dual gluons, 
monopoles and current. 
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' We show that this model exhibits confinement if a condition on the monopole self-energy is 

, satisfied. 
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^ ■ !• INTRODUCTION 

\o . 

In the analysis of lattice gauge theories, it proves very useful to transform between different representations of the 
same theory. In the 70's, Banks, Myerson, and Kogut have shown that U(l) lattice gauge theory can be transformed 
exactly into a representation by a dual photon field and monopoles This photon- monopole representation was 
O ' derived earlier by a different method by Polyakov @, Q . 

One can demonstrate with it that electrostatic charges are confined in 3 dimensions [!, H[ , and that there is a phase 
^3 transition in 4 dimension [H, 0, 0| : in d = 3 the monopoles condense along a string between the charges and create a 
linear confining potential. In dimension 4 there is a critical size of the coupling at which such a condensation sets in. 

The example of U(l) fostered hopes that one could generalize this scheme to non-abelian gauge groups: it led to 
the conjecture that confinement in U(l) and SU(N) result from the same mechanism, and that it has an analogy in 
dual type II superconductors @, lH ■ 

The problem, however, proves to be much harder than in the abelian case. Presently, we only know of three 
non-abelian analogues of U(l) lattice representations: apart from the defining path integral, we have 1. a first-order 
representation, which can be viewed as a lattice version of BF Yang-Mills theory, 2. the spin foam representation 
[9l fiol [TTl H2I ], and 3. the plaquette representation [IH, [lj, [T^]. The first two representations are defined for any 
dimension d > 2, while the latter has been only constructed in 3 dimensions so far. 

It is difficult to go beyond these three representations: both the plaquette and the spin foam representation are 
governed by non-abelian generalizations of constraints that we find in the abelian case. For U(l) these constraints can 
be solved exactly, so that one can move on to obtain the dual Z gauge theory and the photon-monopole representation. 
In the non-abelian case, we do not know at present, how one could solve these constraints exactly, or if there is at all 
a meaningful sense in which representations beyond the known ones exist. 

This article is the third in a series of three papers, where we derive non-abelian analogues of the photon-monopole 
representation for gauge group SU(2). We are not able to do this by exact transformations, and have to rely instead 
on stationary phase app roximations . This leads us to three versions of a gluon-monopole representation of SU(2) 
lattice gauge theory jl6l.[l7j. In the present paper, our starting point is the plaquette representation in 3 dimensions, 
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as it was given by Borisenko, Voloshin and Faber 1 . These authors developed a weak-coupling perturbation theory for 
the plaquette representation, and derived the associated generating functional. It is a path integral over continuous 
and discrete variables that resemble the dual photons and monopoles of U(l). We interpret these variables as a dual 
gluon and monopole field of SU(2). 

We extend the work of Borisenko et al. by including a Wilson loop as a non-trivial background of the perturbation 
theory: we consider the expectation value of a Wilson loop, and treat the latter in such a way that it contributes 
at zeroth order. This is done with the help of the Kirillov trace formula. We then retain only the zeroth order and 
arrive at a gluon-monopole model that has a coupling to a source current J. The structure of the resulting action is 
analogous to the photon-monopole action: it reproduces roughly the tree-level Coulomb interaction, and the coupling 
between monopoles, dual gluons and current is similar to the abelian case. There is an important difference, however: 
it consists in the fact that the dual gluon field is su(2) ~ K 3 -valued (and not M-valued) and that the monopoles couple 
to the length of field vectors. As a result, the gluon-monopole coupling is nonlinear. 

The similarity to U(l) naturally suggests that one could try to generalize Polyakov's method of deriving confinement. 
We propose a way to do so, but it requires an additional heuristic assumption: since we are not able to integrate 
over the dual gluon field as it was done for U(l), we have to assume a renormalization group step that generates a 
monopole self-energy at a scale M below the cutoff scale. With that assumption, the rest follows as in Polyakov's 
paper and one arrives at a non- vanishing string tension for the Wilson loop in the representation j — 1/2. 

The paper is organized as follows: the definition of the plaquette representation is given in sec. [TTJ In sec. IIII1 
we derive the approximative representation by dual gluons and monopoles. The tentative derivation of the area law 
follows in sec. IIVI The results are summarized and discussed in the final section. 

Notation and conventions 

K denotes a 3-dimensional hypercubic lattice of side length L with periodic boundary conditions. The lattice 
constant is a. Depending on the context, we use abstract or index notation to denote oriented cells of n: in the 
abstract notation, vertices, edges, faces and cubes are written as v, e, / and c respectively. In the index notation, we 
write x, (x/i), (x/iu), (xjAVp) etc. Correspondingly, we have two notations for chains. Since the lattice is finite, we 
can identify chains and cochains. As usual, d, d and * designate the boundary, coboundary and Hodge dual operator 
respectively. Forward and backward derivative are defined by 

^fifx = — (fx+afi fx) 7 ^[ifx — ~~ {fx fx—afi) (1) 

where (i is the unit vector in the /z-direction. The lattice Laplacian reads 

A = V M V M . (2) 
For a given unit vector u = fi and a 1-chain J x „, we define 

(u ■ V) -1 J xp .:= ^2 J (xi,...,x'^...,x 3 )f, ■ (3) 

Color indices are denoted by a, b, c, . . . We employ units in which h — c — 1 and a = 1. For some quantities, the 
a-dependence is indicated explicitly. 

II. PLAQUETTE REPRESENTATION OF 3-DIMENSIONAL SU(2) LATTICE YANG-MILLS THEORY 

The basic idea behind the plaquette (or field-strength) formulation is a change of variables from holonomics along 
edges to holonomies around plaquettes. The path integral over the new plaquette variables is constrained, since 
plaquette holonomics have to satisfy Bianchi constraints: the well-known Bianchi identity of abelian lattice gauge 
theory, or generalizations thereof for non-abelian gauge theories. 

There exist different schemes for constructing such plaquette representations (see (13, [lil Il5j ) . In this paper, we will 
use the formulation of Borisenko, Voloshin and Faber for SU(N) lattice gauge theory in 3 dimensions [15J. Although 



1 It is related to earlier plaquette (or field strength) representations by Halpern and Batrouni 1 1. . I lL ITR. 



Figure 1: Eight cubes around an even point: the thick lines indicate a possible choice of connectors. 



the underlying idea is simple, a precise description of this representation is quite involved, as it requires numerous 
conventions on orientations and orderings. For that reason, we will review the essential definitions in this section. 

To simplify the presentation, we will ignore all boundary effects. In the next section, we will restrict ourselves to 
the first term in a weak-coupling expansion, and to that order any boundary-related modifications drop out. The 
proper treatment of boundary conditions is given in [l5[ . 

Edge (or link) variables are denoted by U e , and Wf designates the holonomy around a face /: 

Wf = UUe (4) 

eCdf 

It is assumed that we have chosen some starting point v C /, so that Wf is the product of edge holonomies U e , 
starting at v, and following the orientation of the face /. Whenever we have such (or similar) products of group 
elements, we use left multiplication, i.e. 

n 

Y[U i = U n U n ^ 1 ...U 2 U 1 . (5) 

2 = 

The convention for starting points will be fixed further below. 

In the standard formulation of lattice gauge theory, the expectation value of a Wilson loop C is defined by the path 
integral 



f ( n dC/e 

J \eCK 



(tijWc) = / n dC/ e exp 



-£ ^{Wf + Wf 



4 

/Ck 



Xi (W c ) ■ (6) 



Here, the exponent is given by the Wilson action and 



P=— 2 - (7) 

Wc stands for the holonomy around the Wilson loop C, where, again, a starting point along C is assumed. \j = 
is the character in the representation j. It is understood that for each pair of edge orientations e, e _1 , we integrate 
only over one edge variable U e , with the other one being fixed by 

U e -i =U~K (8) 

The plaquette representation arises from a change of variables: from the edge (or link) variables U e to the Wf's, 
which we call face (or plaquette) variables. This change of variables is achieved in four steps: 

1. Introduce integrations over plaquette variables Wf + delta constraints ^ on them. 

2. Reexpress some of these constraints as non-abelian Bianchi identities. 

3. Impose the maximal axial gauge. 



Figure 2: Cube of k with chosen connector. 



4. Integrate out the edge variables U e . 

The result is a path integral over plaquette variables where of all the constraints on Wf only one type is left, namely, 
the non-abelian Bianchi identities. 

To write down the non-abelian Bianchi identities we need to introduce a number of conventions: consider the set 
of vertices v whose coordinates (x\, X2, x$) are all even and call them the even vertices. Likewise, we call the vertices 
with odd coordinates odd. Clearly, each cube of k contains an even and odd point as corners, and each face of n 
contains either one even or one odd vertex. We will use this to specify our convention for starting points of face 
holonomies Wf : when associating a holonomy to a face we will always take the even or odd point in it as the starting 
point. 

In every cube of k we choose an oriented path consisting of three edges that connect the even and odd point. We call 
this path a connector. We can choose the connectors in such a way that there are altogether four types of connectors 
on the lattice (see Fig. [T]). To each connector we associate a holonomy 

K v , v = U e3 U e2 U ei . (9) 

Let us now specify the non-abelian Bianchi identities. Consider one cube c as in Fig. [2] Take the boundary dc and 
choose an orientation for it. This will induce an orientation for each face / C dc. Call the starting point of the 
connector v and the end point v' . We have three faces of the cube that meet at v: order them fi,fi,f§, starting 
with the face whose orientation agrees with that of the connector, then the one which does not intersect with the 
connector, and then the remaining one (whose orientation is opposite to that of the connector). Likewise, we have 
three faces meeting at v', and we order them /{, fy, J3, according to the same rule as before (see Fig. [2]). 
With these conventions the non-abelian Bianchi identity looks as follows: 

Vc = K-,l W S , Wf, Wf, K p , p W f3 W h W fl =1 (10) 

Wf i denotes the holonomy around the face fi with the starting point at p, and likewise the starting point for Wf; is 
taken to be p' . The reader may convince himself by a drawing that this identity is indeed satisfied. 

In the abelian case, the connectors would cancel and the expression would reduce to the usual abclian Bianchi 
identity. Here, this does not happen. At first sight one might think that this makes a formulation in terms of 
plaquette variables impossible, since the appearance of connectors prevents us from eliminating the edge variables 
completely. It turns out, however, that by the steps 3 and 4 above, the edge variables appearing in (fl"0|) will be 
expressed in terms of plaquette variables, and thus one ends up with a path integral that contains only plaquette 
variables plus constraints on them. 

The expressions for the edge variables depend on how the edge is situated relative to the origin and the even and 
odd points: for example, if x is an even point and x\ > 0, the edge variable U x i is given by (see Fig. [3|) 



U( Xl ,x a ,x 3 ),l = n W(xi,4,0),21 W( XuX2lX > a ) l31 ■ (11) 

x' 2 — X2 — I — 1 

Similar equations hold for other edge variables. To write them down explicitly, we would have to distinguish many 
different cases, although the principle is the same as in equation (fTTj) . In this paper, these details do not matter, since 
we only use the zeroth order of a weak-coupling expansion, and the formulas for the connectors and edge variables 
drop out. 
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Figure 3: Edge variables as products of plaquette variables: in the maximal axial gauge the edge variable U x \ can be expressed 
as a product of plaquette variables. The dots represent even points, and the arrows indicate the order of the product. All edge 
variables on dashed lines equal / due to the maximal axial gauge. 



In a similar fashion the Wilson loop can be expressed as a product of plaquette variables. We take C to be a 
rectangular loop in the 1-2-plane with anti-clockwise orientation, where the lower-left and upper-right corner are 
given by even points (0, 0, 0) and (T, R, 0), T, R > 0. Then, the holonomy around C with starting point x = (0, 0, 0) 
equals 



w c= n n ww^o)^ 



(12) 



xi—T—l X2 — R—0, 



i.e. it is an ordered product of plaquette variables whose plaquettes fill the surface enclosed by the Wilson loop. 

By going through the steps 1 to 4 and using the formulas (HHJ), (fTTj) and (TT2"1) . the original path integral © is 
rewritten as a constrained path integral over plaquette variables Wf. 



(tr 3 W c ) 



(n/^) exp Ef fa ( w > +W 7 1 



Xj (Wc) 



If we change to an index notation, writing (xfj,) instead of e and (xfw) in place of /, the same reads 



(tv 3 W c ) 



I ( n dwv) (ruG4))ex P 

J V XLIV I \ X J 



X f-L<V 



EE - A ^{w^+w-^ v ) 



Xj {Wc) 



(13) 



(14) 



III. REPRESENTATION AS DUAL GLUONS AND MONOPOLES 



After deriving their plaquette representation, Borisenko et al. continue with the construction of a weak-coupling 
(large (3) perturbation theory. The associated generating functional can be considered as a partition function of dual 
gluons and monopoles of SU(2) lattice gauge theory, and generalizes the photon-monopolc partition function of U(l). 

What we add to this scheme is the following: we use the Kirillov trace formula fl9j |. as in paper I, to incorporate 
the Wilson loop in the exponent. With this step we find a generalization of the photon-monopole action that does not 
only contain the dual gluons and monopoles, but also the current of the Wilson loop. This will be crucial in section 
IIVI where we propose a non-abelian generalization of Polyakov's derivation of the area law. 

We will now describe the steps that lead to the description in terms of dual gluons and monopoles: first we will 
rewrite measure, action and Bianchi constraint in terms of Lie algebra elements. To do this for the Bianchi constraint, 
the delta function is expressed as a sum over characters, and the character, by the Kirillov trace formula, as a function 
of a Lie algebra element. We apply the same formula to the trace over the Wilson loop, so that it is parametrized by 
a Lie algebra element. 

In the second step, we apply the Poisson summation formula which trades the discrete sum over representations for 
new degrees of freedom that can be regarded as dual gluons and monopoles. The third and final step is a stationary 
phase approximation. 
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Kirillov trace and Poisson summation formula 

We start by rewriting everything in terms of Lie algebra elements: if we parametrize the plaquette variables as 

W xsxv = j<**°*t* , \ Uxiu ,\ < 2tt, (15) 

the action becomes 

f tr (W XIU , + W-'„) = | cos (|a; x ^|/2) . (16) 

The measure takes the form 

[ <\W - 1 f H 3 ,, sin 2 (|^l/2) ( 

J dW ^- - J {M/2 f (17) 

The Bianchi constraint is written as 

<5(K) = ^(2j + l)xi(K). (18) 

3 

Let v x a a /2 be the Lie algebra element associated to the constraint V x , i.e. 

V x =c lv > a l 2 (19) 

Then the Kirillov trace formula gives 

XjK x> 47rsin(|^|/2) J s2 V ' 

The integral runs over unit vectors n in R 3 , i.e. over the 2-sphere. Next we use the Poisson summation formula to 
replace the sum over representations by an integral and a sum: 

= IV fc >*l/ 2 [ dn e ikn-v x /2 

2 hi 4 ™(M/2) Jsp 

-- I [°° r\rr 2 "V Kl/ 2 f J„ irn-v x /2 + 2m rm 

~ iJ-oo h 4 - sin (Ki/ 2 ) J* 

sin (M/ 2 ) 



R3 



In the last equation we combined the integral over the radius and the unit vector into an integral over R 3 . 
To deal with the Wilson loop, we apply the Kirillov trace formula a second time: by writing 

W c = e iw c°-°/2 (22) 

we obtain 

Xj {Wc) = y + ^l^l{ 2 / dn e 1 ^ ™° / 2 . (23) 

AjV ° J 47rsin(|^ c |/2) J s2 V ' 
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With all this the path integral becomes 



2 I (n/—)/(n dV -)S /«» 



n 



x exp 



B2A0) 

\v x \/2 \ (2j + l)|w |/2 



Bin(|v x |/2)y 47rsin(|w c |/2) 

( 5^/3 cos (K M „|/2) + i(2j + l)n-u; c /2 
2; \m<^ 



+ iip x • w x + 47ri m a 



(24) 

The integration range of uj xfll/ is the ball of radius 27r in R 3 . The scalar fields <p and m are M 3 - and Z- valued 
respectively. Field-independent constants that appear both in the numerator and in Z were dropped. 



Stationary phase approximation 

Let us expand all quantities in (|24| in powers of uj xi1v : the Wilson action yields 

/3cos (|^ Ml/ |/2) = ( 1 - i w 2 ^ + o{cj 3 x ^) 



(25) 



When we expand the non-abelian Bianchi identity (|10p . the contribution of the connectors cancel to lowest order, and 



It follows from equation (TT2"|) that 



IOC 



^ 9 ^xnv^x^v + °{^lpu) 



where S 1 is the minimal surface spanned by the Wilson loop C. From the measure factors we obtain 



\^ X p V \/2 

sin(|w XMIy |/2) 

KJ/2 
sin(|u x |/2) 

sin(|w c |/2) " 



exp 



exp 



1 f\Uxpu\ 



1 

180 
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exp - 



(26) 
(27) 

(28) 
(29) 
(30) 



Each of the terms (|25[) - (|3"0|) contributes with linear, quadratic and higher orders of 

We now apply a stationary phase approximation at u) X nv — and only retain terms up to quadratic order in u) x)Ul . 
We consider the continuum limit where /3 > 1: then, the quadratic term in (|25[> is much larger than the quadratic 
terms coming from ([26l) -(f30" |) . so we keep only the former. This yields the following path integral 2 : 



faWc) = 



n d 3 



x exp 



-82,(0) 

L 

16 



■nix & 



2.7 + 1 

47T 



dn 



sr _p_ 2 it 



(31) 



2 Throughout it is assumed that we apply the same steps within the partition function Z by which we divide. 
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Here, b xiiv stands for the R 3 -valued 2-chain 



b xfiu = (j + 1/2) n S xliU , 



(32) 



i.e. b is proportional to the tensor product of the 2-chain S XI1V (defining the surface S) and the unit vector n in the 
Lie algebra vector space R 3 . We can express it also as 



b xiiv = -u^ (u • V) 1 J xv + u v (u • V) 1 J Xfi 
where u is the unit vector in the xi-direction and the current is defined by 

Jxfi ■= {j + l/2)nC xll . 
Next and integrate over it: 

2.7 + 1 



(33) 



(34) 



(t^c) = |/(ndv)E ^/ 



drc 



x exp 



t \ 2 



2 

E ( ~~q + +47ri|y x |m x 



In this expression, we switched from the 2-chain b xtiv to the 1-chain 

1 



(35) 



(36) 



As in U(l) lattice gauge theory, we can factor off the Coulomb energy by making a change of variables and using the 
identity 3 



(37) 



The final result reads 



2j + l 

47T 



d?i 



x exp 



E ( 75 <^ A< ^ + 4?ri \tp x + A 



, 777 J 7 Q A _1 J a 



(38) 



It should be kept in mind that b xyal depends on the direction of the unit vector n via eq. (|32|) . 

We propose (|38[) as a non-abelian generalization of the photon-monopole representation. The field ip is interpreted 
as a dual gluon field: it has 3 degrees of freedom per point, which agrees with the fact that in 3 dimensions we have 
1 physical degree of freedom per gluon and altogether 3 gluons for SU(2). It mediates the Coulomb interaction 



V. 



jj 



-| E J ^ 1J U = -\ a 9 2 d + V2) 2 E C *pKyC w , 3*0. 



(39) 



The latter agrees roughly with the tree-level result of standard perturbation theory [2(| [2l[ : there one would have 

-\ ag 2 j(j + 1) E C a ^C w . (40) 



T/trcc 

V JJ — ~- 



The discrete variables in arise from the compactness of SU(2) and can be seen as a generalization of monopoles to 
3-dimensional SU(2) lattice gauge theory. 



Eq. H37I I can be proven by starting from the expression e PIJ/V V^6 X „ A 1 e pK \\l and observing that J has no divergence. 
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IV. DERIVATION OF THE AREA LAW 



The close similarity between the photon-monopole representation and (|38[) suggests that it could provide a way to 
generalize the derivation of confinement to 3d SU(2) lattice gauge theory. The main difference consists in the fact that 
the dual gluon field <p is R 3 - and not R-valued and that there appears a modulus in the second term. This renders 
the field theory nonlinear. As a result, we cannot perform a Gaussian integration to write down the analogue of a 
Coulomb gas representation. 

In the abelian case, one can go to the Coulomb gas representation, extract a monopole self-energy and then transform 
back to an effective photon-monopole theory with a lower cutoff scale M. 4 The explicit appearance of the monopole 
self- interaction is important for deriving confinement, since it dampens monopole excitations. 

In the present case, the nonlinearity forbids a simple transition to a Coulomb gas, and we proceed instead as 
follows: we will show that confinement can be derived if the monopole self-energy behaves similarly as for U(l). We 
will assume that after renormalization down to a suitable cutoff scale M, 

• the field (p is replaced by an effective field with a regulated Laplace operator, and 

• a self-energy for the monopoles is generated. 



If this is true, we will get 

fawc) = \ J (n dVx) E 



gj + 1 

47T 



dn 



x exp 



E t; <Px^<Pl ~ 2n 2 (3v ml + 4>ni\<p x + A^Vpfe, 



— T a A -1 J a 

Q XfJ, Xfl 



where the regulated Laplace operator is defined by 

-A" 1 := -A" 1 - f-A 



(41) 



(42) 



and vq is some constant. From here on one can follow Polyakov's arguments to derive the area law Q: let us abbreviate 



r\ x = 4ttA 



fl u Xfl 5 



(43) 



and apply a rescaling and shift on <p 
(tr 3 W c ) = \ 



/(n^)E^/ 



dn e- VlJ 



x exp 



E 



8tt 2 /3 



(fx - Vt^Mpx - Vx) - 2vr 2 /3w m 2 x +im x \(p x 



(44) 



Then, the dilute gas approximation yields 



(tr 3 W c ) 



x exp 



1 

Z 



2j + l 

47T 



dn e 



-Vjj 



E (g^g (V- " Vx)A(fx - r, x ) + 2 e- 2 * 2f3v ° cos |^|) 



Recall that 



J XjJ,V 



= (j + 1/2) nS xll „, 



(45) 



(46) 



4 see the appendix of ref. [IJ and the clarifying remarks in and Q . 
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and observe that the action is invariant under global rotations 

y' a = R\v c , t a ^R a d c xtl . (47) 

The saddle points are determined by the equation 

-A^-^Mf^sinM (48) 

m 

where 

M 2 D = 87r 2 /3e- 2 ' r2/3, ' . (49) 
We replace the regulated Laplacian by the full Laplacian, and plug in eq. (|4"3"| : 

A^ a = 47rV^ - Mf, j— sin |^| (50) 

For simplicity, we assume now that j = 1/2. Consider first the equation in the region above or below the surface S: 

LD a 

Atp a = -M 2 D j— sin | <p\ (51) 

m 

We can find a simple solution for this if we assume that the direction of ip is constant. Then the equation reduces to 
the nonlinear Debye equation for \ip\. 

A\cp\ = -M%sm\(p\, (52) 
We treat this as a quasi 1-dimensional problem and approximate it by the continuum equation 

d 2 M 



-Mlsm\<p\. (53) 



The term 47rV Al 6^ M in ||5D|) is only nonzero at the surface S and gives 

- 4ttV m 6^ = 4tt (4 3 ,o - S X3 - a ) n a . (54) 

This implies that at S the field value ip has to make a jump by Ann. 

If the jump was 27rn, we could construct the solution as in the abelian case for charge q = 1. The fact that it is 
Ann creates a slight (but harmless) complication, and can be treated like the doubly charged loop of U(l) [23j. 

Recall that we have a certain freedom in choosing the particular solution b Xilv . Instead of using the minimal surface, 
so that 

b = nS, (55) 
we could take two surfaces S + and 5_ s.t. dS+ = dS^ — C, and set 

b~ xiiv = (S+ + S-) . (56) 

Now the ip- field has to jump two times by 27rn: once along S-, and a second time along S+. Thus, we can treat the 
situation near each surface similarly as for charge q = 1. 

Imagine that S+ results from "stretching" the minimal surface S to an £3-value W3 > 0, and that S- is the mirror 
image of S+ w.r.t. to the xi-a^-plane. Then, a solution to ([501) i s approximatively given by two domain walls: namely, 

An arctan^-^^ 3 ^ 3 )) - 2tt , (x u x 2 ) & S , x 3 < , 
tp cl (x) w ^ 4n arctan(e- MD ( :E3 ^ 3 )) , (x 1: x 2 ) £ S , x 3 > , (57) 

, otherwise . 
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We plug this solution back into (|45|) and use the trivial saddle point for Z: 



(tVjWc) = 



4-7T 



dn e 



-Vj 



exp 



COS \ (f 



(58) 



The associated action equals AS2dw, where A is the area of the surface S and 6>2dw is the action (or minus the energy) 
of the two domain walls, i.e. 



->2dv 



U1 D 



(59) 



(60) 



(61) 



As in the case of the Coulomb potential, it would be wrong to set j to zero in formula (|46|) . This would yield a 
confining potential when there is actually no Wilson loop. The reason is that we use the Kirillov trace formula in 
conjunction with a stationary phase approximation, and that creates a wrong offset for the j-dependence. The correct 
procedure is to set j — at the beginning of the derivation. Then, the resulting potential is zero, as it should. 



The action does not depend on n, so the integral over n drops out. The result is an area law 

faWc) » 2 exp (-Vjj - aA) 

with a string tension 



V. SUMMARY AND DISCUSSION 



In this paper, we have derived an approximative model for 3-dimensional SU(2) lattice gauge theory. Its degrees of 
freedom can be viewed as a dual gluon and monopole field. We propose it as a generalization of the photon-monopolc 
representation of Polyakov 0, Q and Banks, Myerson and Kogut [l[. 

Our derivation extends an earlier work by Borisenko, Voloshin and Faber [l5j ]. where a weak-coupling perturbation 
theory of the plaquette representation was developed. We only retain the zeroth order of this expansion, but include a 
non-perturbative contribution from a Wilson loop. Thus, our method can be seen as a stationary phase approximation 
around non-trivial field configurations: these are given by monopoles that arise from the compactness of the gauge 
group, and by a worldsheet that is bounded by the Wilson loop. 

The resulting model contains two interaction terms: firstly, a current-current potential that is essentially the 
tree-level Coulomb interaction one would get from a purely perturbative treatment. Secondly, a coupling between 
monopoles, dual gluons and current. This coupling is similar to the photon-monopole coupling of U(l), but nonlinear. 

The analogy with U(l) suggests a possible derivation of confinement. It requires us, however, to make an additional 
assumption on the monopole self-energy. In the abelian case, one can go to the Coulomb gas representation, extract a 
monopole self-energy and transform back to an effective photon-monopole representation. Here, we cannot compute 
the analogue of a Coulomb gas, since the integral over the dual gluon field is non-Gaussian. Instead we assumed that 
a renormalization generates a monopole self-energy at a lower energy scale M. From there on, we can proceed as in 
Polykov's derivation and arrive at a non-vanishing string tension for the Wilson loop. Further investigation has to 
show if our assumption is justified. 

In the two companion papers, we start from two other representations of 3d SU(2) lattice gauge theory and arrive 
at models that are quite similar to the present one. 

Whether the model of this paper is a good or a bad approximation can be tested: it is simpler than the full lattice 
gauge theory and may be easily implemented on a computer. By summing over the monopole variables in expression 
P5[) one can remove the phase factors and translate them into a constraint. The latter can be enforced by a Gaussian 
damping factor. 
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